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Recent Monte Carlo simulation results for the contact values of polydisperse hard-sphere mixtures at a hard planar wall are considered in 
the light of a universality assumption made in approximate theoretical approaches. It is found that the data seem to fulfill the universality 
ansatz reasonably well, thus opening up the possibility of inferring properties of complicated systems from the study of simpler ones. 

+-* 

a . 

g , 1 Introduction 

""O ■ In hard-sphere systems, the general statistical mechanical relation between the thermodynamic properties 
^ \ and the structural properties takes a rather simple form. Since the internal energy in these system reduces to 
q ■ the one of the ideal gas and the pressure equation only involves the contact values of the radial distribution 
l— ~~ '. functions (rdf), a knowledge of such contact values is enough to obtain their equation of state (EOS) and 
| all their thermodynamic properties. However, such a program can not be carried out analytically due to 
J> ■ the lack of exact expressions for these contact values up to the present day. Under these circumstances, the 
^ ! best one can do is to rely on sensible (approximate) proposals based on as sound as possible theoretical 
qq ' results or to rely on computer simulation values. Clearly, the situation is rather more complicated for 
t^. . mixtures than for a single component fluid and in fact for this latter many accurate (albeit empirical) 
equations of state have appeared in the literature from which the contact value may be readily derived. 

A key analytical result is due to Lebowitz [1], who obtained the exact solution of the Percus-Yevick 
(PY) equation of additive hard-sphere mixtures and provided explicit expressions for the contact values of 
' the rdf. Also analytical are the contact values of the Scaled Particle Theory (SPT) [2,3]. Neither the PY 
^ ■ nor the SPT lead to accurate values and so Boublfk [4] (and, independently, Grundke and Henderson [5] 
[ and Lee and Levesque [6]) proposed an interpolation between the PY and the SPT contact values, that we 
■ will refer to as the BGHLL contact values, which leads to the widely used and rather accurate Boublfk- 
O ■ Mansoori-Carnahan-Starling-Leland (BMCSL) EOS [4, 7] for hard-sphere mixtures. Refinements of the 
^ | BGHLL values have been subsequently introduced, among others, by Henderson et al. [8], Matyushov and 
Ladanyi [9] and Barrio and Solana [10] to eliminate some drawbacks of the BMCSL EOS in the so-called 
colloidal limit of binary hard-sphere mixtures. On a different path but also having to do with the colloidal 
l-H ' limit, Viduna and Smith [11] have proposed a method to obtain contact values of the rdf of hard-sphere 
mixtures from a given EOS. In previous work we have made proposals for the contact values of the rdf 



* Corresponding author. Email: malopez@servidor.unam.mx 
t Email: santos@unex.es 
^ Email: andres@unex.es 

Molecular Physics 

ISSN 0026-8976 print/ ISSN 1362-3028 online ©2005 Taylor & Francis Ltd 
http:/ / www.tandf.co.uk/journals 
DOI: 10.1080/002689700xxxxxxxxxxxx 



2008 



21:6 Molecular Physics Liblice06'rev 



2 Lopez de Haro, Yuste and Santos 

valid for mixtures with an arbitrary number of components and in arbitrary dimensionality [12, 13] and 
for a hard-sphere polydisperse fluid [14], that require as the only input the EOS of the one-component 
fluid. Apart from satisfying known consistency conditions, they are sufficiently general and flexible to 
accommodate any given EOS for the one-component fluid. As far as computer simulation results are 
concerned, contact values of the rdf of hard-sphere systems have been reported by Lee and Levesque [6], 
Barosova et al. [15], Lue and Woodcock [16], Cao et al. [17], Henderson et al. [18], Buzzacchi et al. [19] 
and Malijevsky [20]. In particular, these last two references study hard-sphere mixtures in the presence of 
a hard planar wall. 

It is interesting to point out that in the case of multicomponent mixtures of hard spheres and in the 
polydisperse hard-sphere fluid, the contact values which follow from the solution of the PY equation [1], 
those of the SPT approximation [2,3], those of the BGHLL interpolation [4-6] and our own prescriptions 
[12-14] exhibit a feature that one might catalogue as a "universal" behavior because, once the packing 
fraction is fixed, the expressions for the contact values of the rdf for all pairs of like and unlike species 
depend on the diameters of both species and on the size distribution only through a single dimensionless 
parameter, irrespective of the number of components in the mixture. 

In Ref. [14] we have presented a comparison of the different theoretical proposals for the contact values 
of the rdf and the ensuing EOS stemming from them with the simulation results. The aim of the present 
paper is to assess, apart from the accuracy of such proposals, whether the universality feature alluded to 
above is indeed present in the simulation data in the case of mixtures in the presence of a hard wall. This 
represents an extreme case and therefore a proper test ground for our approach. 

The paper is organized as follows. In order to make the paper self-contained, in Sec. |2]we rederive our 
most recent proposal [14] for the contact values of the rdf (labelled as e3 for the reasons explained below) 
using some known consistency conditions and two different routes to compute the compressibility factor 
of a polydisperse hard-sphere system in the presence of a hard planar wall. We also point out in this 
section that the other proposals sharing the universality feature, namely the PY, SPT, BGHLL and our 
two previous proposals [12, 13] may be cast in the same form as our e3 proposal, but that only the SPT and 
the e3 proposals are consistent in the sense that they lead to the same compressibility factors with the two 
different routes. Section 01 deals with the comparison between the various contact values and simulation 
results, examining both the accuracy of the theories as well as whether the universality ansatz is confirmed 
by the simulation data. We close the paper in Sec^Jwith further discussion and some concluding remarks. 



2 Contact values of the radial distribution functions 

Consider a polydisperse hard-sphere mixture with a given size distribution f(o~) (either continuous or 
discrete) at a given packing fraction rj = ?p/^3, where p is the (total) number density and 

/■oo 

t i n = (a n )= / daa n f(a) (1) 
J o 

is the n-th moment of the size distribution. Let g(a, a') denote the contact value of the pair correlation 
function of particles of diameters a and a'. This function enters into the virial expression of the EOS 
as [21] 

Z ee = 1 + O- r da da' f{a)f{a>) ( Z+^L) ' g{a , a') 

pk B T p 3 J J V 2 / 

= l + ^((a + afg(a,a')), (2) 

where Z is the compressibility factor, p is the pressure, ks is the Boltzmann constant and T is the absolute 
temperature. Assume further that the polydisperse hard-sphere mixture may find itself in the presence 
of a hard wall. Since a hard wall can be seen as a sphere of infinite diameter, the contact value of the 
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correlation function g w (cr) of a sphere of diameter a with the wall is obtained from g(a,a') as 

g w {a) = lim g{a,a'). (3) 

Note that g w (cr) = p w (cr) //Obulk (c) provides the ratio between the density of particles of size a adjacent to 
the wall, p w (a), and the density of those particles far away from the wall, PbulkC ")- There is a sum rule 
connecting the pressure and the above contact values [22], which provides an alternative route to the EOS, 
namely 

/>oo 

da f(a)g w {a) = (g w (a)), (4) 

where the subscript V0 111 Z yj tlclS been used to emphasize that Eq. (jlj) represents a route alternative to the 
virial one, Eq. (J2J), to get the EOS of the hard-sphere polydisperse fluid. Our problem is then to compute 
g(a,a') and the associated g w (a) for the polydisperse hard-sphere mixture in the presence of a hard wall, 
so that the condition Z = Z w is satisfied. 
We consider a class of approximations of the type [13, 14] 

g(a,a')=G(z(a,a')), (5) 

where 

, _ 2aa' [i 2 , , 

z(a,a) = — (6) 

cr + cr'/i 3 

is a dimensionless parameter. Therefore, at a given packing fraction r/, we are assuming that all the 
dependence of g(a,a') on cr, a' and on the details of the size distribution /(cr) is through the single 
parameter z{a, a'). This implies that if two different pairs (<ta, &'a) an d (c\B) &'b) i n two different mixtures A 
and B (at the same packing fraction) have the same value of the parameter z, i.e., za(cta, cr' A ) = zb^ob^'b)-, 
then they also have the same contact value of the rdf, i.e., gA^A,^) = gB(aB,cr' B ). The parameter 
z^ 1 (a,a') = (cr -1 + cr' -1 )/(2/i2/ ^3) can be interpreted as the arithmetic mean curvature, in appropriate 
units, of spheres cr and a' [13]. 

Notice that Eq. © implies in particular that g w {a) = G(z w (a)), where z w (a) = 2cr/i 2 //i 3 . Once one 
accepts the "universality" ansatz ©, the remaining problem lies in determining the form of the function 
G(z). This may be achieved by considering some consistency conditions. Note that in the one-component 
limit, i.e., f(a) = 5(a — do), one has z = 1, so that [12,13] 

G(z = 1) = g p , (7) 

where g p is the contact value of the radial distribution function of the one-component fluid at the same 
packing fraction 77 as the packing fraction of the mixture. Next, the case of a mixture in which one of the 
species is made of point particles, i.e., a — > 0, leads to [12-14] 

G(z = 0) = — !— = G . (8) 
1 — r] 

We now want consistency between both routes to the EOS for any distribution /(cr). To this end, we 
assume that z = is a regular point, take into account condition © and expand G(z) in a power series 
in z: 

CO 

G{z) =G + Y,GnZ n . (9) 

n=l 
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Using the ansatz Q and Eq. Q in Eq. (j2j) one gets 

oo n 

Z = l + r,J2 2 n - 1 G n J^ T (a n a' n (a + a'f~ n ) 

n=0 ^3 

oo n 

= Go + 3^G + V £ 2 n - l G n J^(a n a' n (a + af~ n ), (10) 

^ ti ^3 + 

where in the last step we have taken into account that 

1 + rj—^-((a + a') 3 ) = 1 + n— (/i 3 + 3/ii/x 2 ) 

2/U 3 /i3 

= G + 3^G . (11) 

^3 



Analogously, Eq. @ becomes 



Z w = Go+} 2 n G n —\x n . (12) 
i M 3 

n=l 

Notice that if the series © is truncated after a given order n > 3, Z„, is given by the first n moments of the 
size distribution only. On the other hand, Z still involves an infinite number of moments if the truncation 
is made after n > 4 due to the presence of terms like (a 4 a' 4 /(a + a')) in Eq. (|TU|) . Therefore, if we want 
the consistency condition Z = Z w to be satisfied for any polydisperse mixture, either the infinite series @ 
needs to be considered or it must be truncated after n = 3. The latter is of course the simplest possibility 
and thus we consider the approximation 

G(z) =G Q + G lZ + G 2 z 2 + G 3 z 3 . (13) 

As a consequence, Z and Z w depend functionally on f(cr) only through the first three moments (which is 
in the spirit of Rosenfeld's Fundamental Measure Theory [23]). 
Using the approximation (|T3|) in Eqs. (fTT)|) and (|T2|) we are led to 



Z = G + n 



3 

^ (3G + 2Gi) + 2^| (Gi + 2G 2 + 2G 3 ) 
^3 M 3 



(14) 



Z w = G + 2^G X + 44 (G 2 + 2G 3 ) . (15) 
M3 Ml 

Thus far, the dependence of both Z and Z w on the moments of f(cr) is explicit and we only lack the 
packing- fraction dependence of Gi, G 2 and G3. From Eqs. (|14|) and ()15|) it follows that the difference 
between Z and Z w is given by 

3 

Z - Z w = ^ [37/Go - 2(1 - n)Gx] + [r?Gi - 2(1 - r?)G 2 - 2(2 - r ? )G 3 ] . (16) 
m M3 

Therefore, Z = Z w for any dispersity provided that 

Gi = (17) 

2(l-r?) 2 ' 
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G2 = —p- — 73 ~ 1 -^3, (18) 

4(1-7/)^ 1-?? 

where use has been made of the definition of Go, Eq. (|SJ). To close the problem, we use the equal size limit 
given in Eq. Q and after a little algebra we are led to 

G 2 = {2-ri)g p - 2 -±^, (19) 
(I-77) 



;i - v) (gf T - sp) 



(20) 



where 



SPT 1 - 77/2 + r? 2 /4 
h (1 _ r? )3 



(21) 



is the contact value of the radial distribution function for a one-component fluid in the SPT. This completes 
our derivation of the e3 approximation leading to the two following main results for the contact values [24] : 



, 1 377 H2 era 1 
g(a,a ) = h — : — 7 + 4 



I-?? (1 - r\) A*3 a + c' 
+8 U -,;,(,r-.%) / " 2 ™' 



2 + ?7 2 /4' 



H2 aa 
Hza + a' 



/j, 3 a + a' 



(22) 



9w(v) 



(1-1?) w 

SPT _ \ I ^2 



1 — 77 

+8(1 - t?) (<7 P 



2 + t? 2 /4 



(i-7?r 



^2 
— < 

/*3 



3p) l —° 



(23) 



The label e3 is meant to indicate that (i) the resulting contact values are an extension of the one-component 
contact value g p and that (ii) G(z) is a cubic polynomial in z. As mentioned earlier, all the theoretical 
proposals that also comply with the universality ansatz namely the PY, SPT, BGHLL and our two 
former proposals [12, 13] for the contact values of the rdf , may be written in the form of Eq. ()13|) with 
G3 = 0, but only the SPT values also yield Z = Z w for any dispersity (see Table I in Ref. [14] for 
details). Note further that the practical application of Eqs. (|22[) and (|23|) needs only the specification of 
the size distribution f(a) and the choice of an approximate expression for g p . For the latter, we will use 
the Carnahan-Starling EOS [25], namely 



cs _ 1-77/2 
y p (! _ r? )3 



(24) 



and use the notation eCS3 to label the approximation. As for the size distribution, we will consider three 
cases: 
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Tabic 1. Parameters of the size distributions for 
the examined mixtures. 



Type Parameters n 



Top- hat 




c = 0.2 




0.2 


Top-hat 




c = 0.2 




0.4 


Top-hat 




c = 0.7 




0.4 


Schulz 




g = 5 




0.2 


Schulz 




q = 5 




0.4 


Bidisperse 


0-2/0-1 


= 3, £'2 = 


= 0.02 


0.206 


Bidisperse 


0-2/0-1 


= 3, X2 - 


= 0.04 


0.207 


Bidisperse 


02/01 


= 3,22 = 


= 0.06 


0.208 


Bidisperse 


0-2/0-1 = 


= 3, x 2 = 


0.0193 


0.404 


Bidisperse 


0-2/0-1 = 


= 3, X2 = 


0.0358 


0.407 


Bidisperse 


0-2/0-1 = 


= 3, X2 = 


0.0621 


0.401 



(i) The top-hat distribution of sizes given by 




I 0, otherwise, 
(ii) The Schulz distribution of the form 

q+lfq+1 V ( q+1 \ 

f(u) = — — a exp a 26 

qlM \ Hi J \ Hi J 

(hi) The case of a bidisperse mixture, namely 

f(a) = xi5(a - cti) + X2d(a - ct 2 ), x\ = \ — xi. (27) 

This choice of size distributions may seem to be to some extent arbitrary (one could for instance have 
also included a log-normal distribution). It has been mainly motivated by our desire to compare with the 
(to our knowledge) available simulation data for polydisperse hard-sphere mixtures in the presence of a 
hard planar wall. Moreover, those simulations have been computed for common packing fractions in the 
polydisperse systems (I25|) - (|27j) . In Tabled we present the values of the parameters corresponding to the 
polydisperse mixtures that have been recently studied, via Monte Carlo (MC) simulations by Buzzacchi 
et al. [19] and the bidisperse mixtures studied also using MC simulations by Malijevsky [20]. We are now 
in a position to assess the merits and limitations of our proposal. 



3 Comparison with simulation results 

In Figs, n and El we display the comparison between the results of various approximate theories for the 
contact values of the wall-particle correlation functions and those obtained from computer simulations for 
the mixtures given in Table ^ [26] . In view of the fact that our main concern is to assess the universality 
ansatz, we have chosen to represent the difference g w {u) — 5^ GHLL as a function of z w (a)/2 = 0^2! The 
figures suggest that, although the simulation data that we have examined are only a few, the universality 
ansatz seems to be followed by them to a large extent, thus providing an a posteriori support to the theo- 
retical approaches that have this feature. In particular, the data corresponding to the several polydisperse 
and bidisperse mixtures with 0*112/ Hz ^5 1-5 overlap reasonably well. In addition, the isolated points cor- 
responding to the data for the big spheres in the bidisperse mixtures are consistent with the trend shown 
by the points with 0^2! Hz ^5 1-5- Of course more simulations for other values of 0H2I Hz, especially in the 
region o~Hil Hz ^ 1-5 or including ternary systems with higher diameter ratios (which would offer a more 
stringent test at higher values of z), would be welcome to further confirm this assertion. 

Regarding the theoretical approaches, as clearly seen in the figures and already mentioned in Ref. [14], 
the overall trend is captured best by the eCS3 approach, while the PY approximation values are even 
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Figure 1. Plot of the difference of contact values g w — <?2 GHLL as a function of 07^2/^3 for the different polydisperse mixtures at the 
fixed packing fraction r\ = 0.2. The symbols are MC simulations: top-hat distribution with c = 0.2 (Ref. [19]) (squares), Schulz 

distribution with q = 5 (Ref. [19]) (triangles), bidisperse mixtures (Ref. [20]) (filled circles). The lines are PY( ), SPT ( — ■ — ■) and 

eCS3 (— ). 




Figure 2. Plot of the difference of contact values g w — 3§ GHLL as a function of 07^2/^3 for the different polydisperse mixtures at the 

fixed packing fraction n = 0.4. The symbols are MC simulations: top-hat distribution with c = 0.2 (Ref. [19]) (squares), top-hat 
distribution with c = 0.7 (Ref. [19]) (circles), Schulz distribution with q = 5 (Ref. [19]) (triangles), bidisperse mixtures (Ref. [20]) (filled 

circles). The lines are PY( ), SPT ( ) and eCS3 (— ). 

qualitatively at odds with the simulation data. The SPT overestimates the contact values in the region 
o'A i 2// u 3 ^ 1 but it becomes the second best approximation for larger values of a^/nz- In the latter region, 
all the theories underestimate the simulation data, the eCS3 predictions being the most accurate, especially 
for the smallest packing fraction. 



4 Concluding remarks 



In this paper we have examined the universality assumption that is present in many theoretical derivations 
by which, once the packing fraction is fixed, for all pairs of like and unlike spheres in a polydisperse hard- 
sphere mixture with an arbitrary size distribution and in the presence of a hard wall, the dependence of 
the contact values of the particle-particle correlation function, g(cr,a'), and of the wall-particle correlation 
function, g w (a), on the diameters and on the composition is only through a single dimensionless parameter 
and holds for an arbitrary number of components. This was done by comparison with available MC 
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simulation results for g w (a) because, since z w (a) > z(a, a'), these contact values represent a more stringent 
test for the universality ansatz than the values of g(a, a'). While our analysis is limited due to the few data 
that are at hand, the results suggest that indeed the simulation data seem to comply reasonably well with 
the ansatz, the results corresponding to the three different bidisperse mixtures virtually falling on top of 
the polydisperse ones for common values of z w . The results also indicate that our eCS3 approximation does 
a rather reasonable job. Although it underestimates the simulation data for high values of the parameter 
z w , it is still better than the other approximations sharing the same universality property. 

A noteworthy aspect of the comparison between the simulation data and the theoretical approximations 
is that those proposals that fulfill the condition Z = Z w , namely the SPT and eCS3, are the ones that 
show the best performance for high z w . Since, as shown by Eq. (|2T))) and discussed in Ref. [14], our e3 
approach becomes identical to the SPT one when the choice g p = g p PT instead of g p = <?p s is made, we 
can interpret it as a versatile and flexible generalization of SPT. We are fully aware that, apart from the 
consistency conditions that we have used, there exist extra ones (see for instance Ref. [27]) that one might 
use as well within our approach. Assuming that the ansatz © still holds, these conditions are related to 
the derivatives of G with respect to z, namely 



8G{z) 



dz 



3// 



2(1-7/)' 



(28) 



d 2 G(z) 



dz 2 



3?/ 



z=0 



gpy - -^9 



(29) 



d 3 G(i 



dz 3 



0.. 



(30) 



where g p Y = (1 + r//2)/(l — r]) 2 is the contact value of the one-component hard-sphere fluid in the PY 
approximation. The question immediately arises as to whether the fulfillment of these extra conditions 
might influence the results we have presented in this paper. Interestingly enough, as shown by Eq. (|17j). 
condition (|28|) is already satisfied by our e3 approximation without having to be imposed. On the other 
hand, condition ()3(Jj) implies G3 = in the e3 scheme and thus it is only satisfied if g p = flp PT , in which 
case we recover the SPT. Condition (|29|) is not fulfilled either by the SPT or by the e3 approximation 
(except for a particular expression of g p which is otherwise not very accurate). Thus, fulfilling the extra 
conditions (j29j) and (|3()|) with a free g p requires either considering a higher order polynomial in z (in which 
case the consistency condition Z = Z w cannot be satisfied for arbitrary mixtures, as discussed before) or 
not using the universality ansatz at all. In the first case, we have checked that a quartic or even a quintic 
polynomial does not improve matters, whereas giving up the universality assumption increases significantly 
the number of parameters to be determined and seems not to be adequate in view of the behavior observed 
in the simulation data. Therefore, e3 appears to be a very reasonable compromise between simplicity and 
accuracy, with the added bonus of being versatile to accommodate any choice for g p . 

Finally, one should point out that the fact that the simulation results give support to the validity of the 
universality assumption, opens up the possibility of gaining information of rather complicated polydisperse 
mixtures from the knowledge of simpler systems using an approximation like our e3 approximation. 

We want to thank Dr. Alexandr Malijevsky for making his simulation results available to us prior to 
publication and Dr. David Reguera for bringing Ref. [27] to our attention. M.L.H. acknowledges the finan- 
cial support of DGAPA-UNAM under project IN-110406. The research of A.S. and S.B.Y. has also been 
supported by the Ministerio de Education y Ciencia (Spain) through grant No. FIS2004-01399 (partially 
financed by FEDER funds). 
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